Evidences of Bolgiano scaling in 3D Rayleigh-Benard convection 
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We present new results from high-resolution high-statistics direct numerical simulations of a tri- 
dimensional convective cell. We test the fundamental physical picture of the presence of both a 
Bolgiano-like and a Kolmogorov-like regime. We find that the dimensional predictions for these two 
distinct regimes (characterized respectively by an active and passive role of the temperature field) 
are consistent with our measurements. 



The dimensional theory of homogeneous and isotropic 
turbulence has been definitely settled long ago by the 
work of A. N. Kolmogorov j^]. On the other hand, it 
is still missing a clear theoretical picture for the strong 
fluctuations in the energy dissipation field that lead to in- 
termittency effects (i.e. non gaussian behaviour of prob- 
ability distribution functions). Phenomenological theo- 
ries have been proposed but no systematic theory for 
computing experimentally measured numbers has been 
successful so far. The situation of "non-ideal" turbu- 
lence is even more controversial, already at the level of 
dimensional expectations. A typical realization, the one 
we will address in this paper, is the tri-dimensional (3D) 
Rayleigh-Benard cell, described, in the Boussinesq ap- 
proximation by the following set of equations: 



dtv -\- (v ■ \7) V = — Vp -f 
dtT +{vW)T = xV^T 



agTz 



(1) 
(2) 



with isothermal boundary conditions on the upper and 
lower planes of a cell of height H: T{z = 0) = -\-AT/2 
and T{z = H) = -AT/2. As usual, vlx,t) is the 
velocity field and T{x, t) the temperature field. Kine- 
matic viscosity and thermal conductivity are respectively 
v and X, while the thermal expansion coefficient is a 
and gravity acceleration is g. In the following we will 
mainly focus on the longitudinal structure functions of v 
and T: Sp{r) = { [{v{x + r) - v{x)) ■ r]^ ) and Tp{r) = 
{[Tix + r)~Tix)r). 

In this work we present some tests of the predictions for 
the structure functions defined above that can be derived 
in the scenario proposed years ago by Bolgiano to 
describe convective turbulence. 

Despite much research on the subject, sound evidence 
of the validity of the Bolgiano scenario and the recov- 
ery of Kolmogorov scaling at small scales is still missing. 
Good quality confirmation of the dimensional Bolgiano 
scenario was recently shown in a two dimensional nu- 
merical simulations p^ . This result cannot be directly 
related to the 3D case because of the strong differences 
in the properties of the velocity field in 2D. Statistical 
properties of the velocity field were recently used with 



even simpler models (shell models for turbulence) p2[ . 
This approach is even less probing since these models 
were built precisely in order to implement Bolgiano scal- 
ing. Finally, a recent experiment even questioned the 
behaviour the the Bolgiano length inside a convective 
cell ll^. Of course from an experimental point of view 
it can be difficult, if not impossible at all, to measure 
all relevant quantities and it can be even harder to have 
access to information at several (not just a few) posi- 
tions inside the volume. The results presented in this 
letter show good consistency with the idea of a Bolgiano 
regime (i.e. a range of length scales where temperature 
driven buoyancy effects are dominant) but, because of 
the still limited resolution of "state of the art" numerical 
simulations, we will have to resort to somehow indirect 
tests. 

This letter is organized as follows: a brief review of 
phenomenological expectations, details of our numerical 
simulations, data analysis and then concluding remarks. 

Starting from equations (|^) and (||), if one uses di- 
mensional analysis and assumes homogeneous scaling for 
velocity and temperature differences (inside the inertial 
range), one ends up with two distinct scaling regimes. At 
small scales (Kolmogorov-like scenario), r <C Lb, 

Sv{r) - £1/3^1/3 

ST{r) - ivi/2e-i/6^i/3 (4) 

while at large scales (Bolgiano-like scenario), r ^ Lb, 

5v{r) ^ {agf'"^ N^/^^'"^ (5) 
5T{r) ^ (agy^'^ N^/\^/^ (6) 

The Bolgiano length. Lb, is an estimate of the distance 
at which the dissipative and buoyancy terms on the right 
hand side of equation (^ balance, valid under the as- 
sumption of a scaling behaviour. In the following we will 
use the z dependent version of Lb{z) introduced in Q| 
that uses averages of e and N defined at a given height z, 
e{z) = ;./2(E„- {d^v. f)^ and iV(z) = x/2(E. (^^Tf).■■ 



Lb{z) 



z{zfl^ 



N{zY/^{agf'^ 



(7) 
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FIG. 1. Log-log plot of the structure function 5'3(r), de- 
fined in the text, measured close to the end plates and 
at the center of the cell (x). The horizontal scale is in grid 
points, while the vertical scale is in arbitrary units. 
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FIG. 2. Plot of the z-dependent Bolgiano length as a func- 
tion of the distance from the bottom wall, z, as given by the 
usual definition of eq. |^ {Lb{z), +) and as given by the pro- 
cedure described in the text, following eq. H {a ■ Lb{z), x). 
All dimensions are in grid points and a = 3.1. 



Our analysis has been performed on data coming from 
Direct Numerical Simulations (DNS) employing a by 
now standard Lattice Boltzmann scheme on a mas- 
sively parallel computer ||^. The resolution of the nu- 
merical simulation was 240^ and the Rayleigh number 
[Ra = agATH^/{vx)) was approximately 3.5 • 10^. The 
Prandtl number was equal to unity. We performed a 
stationary simulation extending over approximately 500 
recirculation times and stored nearly 400 independent 
configurations with full information on all velocity com- 
ponents and the temperature field. Boundary conditions 
were periodic in the x and y directions (in order to main- 
tain homogeneity on horizontal planes) and isothermal at 
the top and bottom planes of the cell {z = and z = H) . 

Free slip boundary conditions (i.e. satisfying only 
the incompenetrability constraint) were used on the 
top/bottom planes for the velocity field. This choice was 
made in an attempt to reduce the effects of a viscous 
boundary layer close to the horizontal walls. Indeed, as 
will be clear from the following, thermal effects are domi- 
nant near the isothermal walls, so using no-slip boundary 
conditions might produce effects on the velocity statistics 
interfering with the ones coming from pure buoyancy. 
More details on this point will be given later on. 

The most direct way to test the dimensional validity 
of the Bolgiano picture would be to measure structure 
functions and to check whether they scale with exponents 
close (apart from intermittency corrections) to the ones 
predicted by the sets of equations (^]4|) and (|[]|). 

Unfortunately this cannot be done directly because of 
the limited resolution of our DNS. 

Fig. |l| substantiates this comment, by plotting 
the velocity structure function of order 3, 5*3(2;, r) = 
( {v{x + r,z) — v{x, z))^ ) at two different position: one 
very close to the wall (z = 10), and the other at the 
center of the cell (z = 120 = H/2). 



As it can be seen from Fig. no evident scaling range 
can be detected, even if a steepening of S'3(r) is clearly 
detected as one comes close to the walls. 

Given this state of affair, in order to test the consis- 
tency with the Bolgiano picture we performed two dis- 
tinct, even if somehow less direct, tests. 
The first test consists in checking that the Bolgiano 
length actually keeps track of the scales at which the 
buoyancy term balances the dissipative term in equation 
®. ■ 

First, we measure Lb{z) in terms of e{z) and N(z) 
and plot our results in Figure From a first look at 
the behaviour of Lb{z) we learn that Bolgiano effects 
should be measurable, if at all, near to the isothermal 
walls (where Lb{z) is of the order of 10^). 

Close to the center of the cell. Lb is of the order of lO'^ 
and Kolmogorov like behaviour is expected to be mea- 
surable at almost all scales. 

We then measure, directly and independently from 
the previous quantity, the scale at which dissipation and 
buoyancy effects balance, i.e. we look for the scale Lb{z) 
for which: 

e{z) ^ ag{5vCLBW{LB)), (8) 

We provisionally consider Lb at modified definition of 
the Bolgiano length. In Fig. ||we plot Lb{z) and aLB{z) 
with a constant which was tuned to be a '--^ 3.1. 

As it can be seen, the two definitions yield the same 
behaviour apart from the multiplicative factor, a. The 
reason for the multiplicative factor (of order unity) is 
due to the fact that the two definitions are dimensional 
estimates so they can miss a numerical prefactor. Con- 
sidering this, the fact that the two definition behave in 
the same way after rescaling has to be considered, in our 
opinion, an excellent agreement. 

Here we want to underline two points which, we be- 
lieve, add relevance to this finding. First of all the 
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two definitions are of course linked but definitely differ- 
ent. The "traditional" definition of Lb (see also [Q), 
as from definition (0), comes from supposing that the 
two scaling laws in equations (||) and (||) merge at Lb, 
hence comes from solving the equation e{zY^^LB{zY/'^ — 

{agf^^Nizf/^LBizf^- The second definition, Lb{z) 
comes instead from a direct measurement of the strength 
of the dissipative and forcing term in equation (Q) . 

The second important point consist on the fact that 
our cell is not homogeneous; this adds strength to the 
equality between the two different definitions of Lb{z). 

As a consequence of this test we can claim that the Bol- 
giano scenario and the expected power law behaviour are 
consistent with the value measured for the terms appear- 
ing on the right hand side of equation (|l|). Furthermore 
we fully confirm (with higher accuracy) our former re- 
sults for Lb . We like to underline that extracting the 
behaviour of the Bolgiano length using measured quan- 
tities that are not those appearing in its definition could 
introduce large errors. For example in [ p^ a Bolgiano 
length was extrapolated as the scale at which there is 
a change of slope in a particular structure function. In 
using such a procedure one is dominated by strong finite 
size effects present on the structure functions. 

We now proceed to our second test. Under the hy- 
pothesis of validity of Bolgiano scaling, and if enough 
resolution were available, one would expect to see power 
law behaviour in the inertial range in Fig. |l] with two 
distinct slopes 1 and 9/5, corresponding respectively to 
Kolmogorov and Bolgiano scaling. With available reso- 
lution, we are not able to detect a clear power law be- 
haviour from Fig. |^, although we see a clear steepening 
of the structure function as we come close to the wall. 
What we are going to do in the following is to try to 
quantify as well as possible this change of slope. 

We adopt the following procedure. We focus on the 
plot of SQ{r) versus 5*3(7-) and apply Extended Self Simi- 
larity (ESS, see 0) in order to detect a trustable plateau 
in the local slopes. We found this plateau to corre- 
spond roughly to distances in the interval X" = [25, 40] 
for the velocity structure function and 1^ = [15, 30] for 
the structure functions of the temperature field. We 
then define two other interval slightly shifted to the left 
II = [20, 35], 2^ = [10, 25] and to the right II = [30, 45], 
— [20,40]. These interval were of course shifted by 
a reasonable amount, i.e. were still possible highest or 
lowest estimate for the same plateau. 

We finally perform a power law fit to extract a scaling 
exponent on the structure functions for the velocity and 
for the temperature in the three intervals defined before. 
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FIG. 3. Behaviour of the exponents C,3,{z), Xai.^) s-s a func- 
tion of z. 

We define the scaling exponents for the structure func- 
tions of interest as follows: 

{5v{z,rY) ^M''^ (9) 
{5T{z,rY) r^r^^{z) (10) 

{{5v{z,r)5T{z,rf)^' ) ^ rP^'^'^ (11) 

Our fits provide a central value for the exponents 
and two, respectively, higher and lower estimates (cor- 
responding to the shifted intervals). The procedure 
adopted in estimating the errors reflects the fact that 
the largest source of systematic error is connected with 
the choice of the fitting ranges and not with statistical 
accuracy. 

In the Kolmogorov regime we expect the exponents of 
equations dgfjll]) to take the following dimensional val- 
ues: Cj, = p/3, Xp — p/3 and pp — p/3. In the Bolgiano 
regime, on the other hand, we expect the following di- 
mensional values: Cp — 3/5p, Xp = p/5 and pp = p/3. 

In Figure || we plot the behaviour of the measured cen- 
tral exponents and their errors as a function of z for C3{z) 
and Xii^)- The behaviour is qualitatively and quantita- 
tively consistent with the expected scaling exponents: we 
observe a smooth transition from a Bolgiano dominated 
regime (near to the wall, small z) to a Kolmogorov regime 
(approaching the center of the cell). 

In Figure ^ we have not shown the behaviour of P3{z) 
as it is consistent (within error bars) with the constant 
value 1 and plotting it on the same figure would have 
made it unreadable. 

Here we, once again, like to underline the consistency 
of our findings. From the theoretical picture we expect 
to see the crossover to the Kolmogorov scaling when 
Lb{z) < H. This is indeed the case as Lb{z) < H for 
z < 40 and indeed we see Csi^) ~ 1 in the same range of 
z. 

Another interesting question concerns the "real" statis- 
tics of the Bolgiano regime. Very recently interest is 
growing along this line of research because of the desire to 
understand the differences of the statistical properties of 



3 



an active with respect to a passive scalar [|ri|,Q . In this 
work we focused on the gross features, i.e. on the dimen- 
sional behaviour. If we try to look at intermittency, by 
means of ESS we find a strong increase of intermittency 
for the velocity field approaching the isothermal walls. 

Unfortunately we believe that with our simulation we 
are not in a position to make any definite statements 
about intermittency in the Bolgiano dominated regime. 

Indeed a study of intermittency in the Bolgiano regime 
would involves positions nearby the isothermal walls (as 
only there the Bolgiano length scale is small enough to 
have an inertial range largely dominated by buoyancy 
effects). Recently it was found that intermittency in- 
creases in the velocity structure functions inside a vis- 
cous boundary layer . In an attempt to reduce the vis- 
cous boundary layer we decided to apply free-slip velocity 
boundary conditions to the isothermal walls free slip for 
the velocity. Unfortunately our idea cannot be successful 
because of two reasons. First, the free-slip boundary con- 
ditions does not completely suppress the boundary layer. 
A boundary layer thickness can be defined through an ex- 
trapolation of {vz{x)'^), and in our case turned out to be 
of the order of roughly 15 grid spacings. Second, it was 
recently realized that mechanical turbulence in a free-slip 
channel also presents an enhancement of intermittency 
near to the walls 

A procedure to disentangle buoyancy and planar effect 
is clearly needed to make any definite statement on inter- 
mittency in the Bolgiano regime. Otherwise one cannot 
decide whether a change of intermittency is related to 
Bolgiano dynamics or to the increase occurring nearby 
boundary layers. 

In order to clarify this point it would be important to 
perform a simulation of a Rayleigh-Benard-like system 
with periodic boundary condition in all directions (i.e. 
an homogeneous Rayleigh-Benard cell). We suggest to 
extend to 3D the study made in 2 dimensions in |l l| . 

In order to do that one could write temperature field 
as the sum of a linear profile plus a fluctuating part, 
T(x, y, z) = T„„(z) + T'(x, y, z), with T„„ = AT/2 ■ (1 - 
2z/H) and obtain: 



Choosing the parameters in order to have a Bolgiano 
length as small as possible, one would benefit of a wide 
range of scales were to study the buoyancy dominated 
flow. Further advantage of homogeneity would be the 
natural increase of statistics and also the applicabil- 
ity of tools like SO{?>) decomposition, to disentangle 
anisotropic terms |l^JT^. A study of this kind is in 
progress. 

Concluding we have performed a number of basic tests 
in order to validate the scenario of Bolgiano-Kolmogorov 



scaling in a convective cell, within the limitations but 
also advantage of z non homogeneity of our cell. We were 
able to confirm the transition between the two expected 
scenarios. 

We acknowledge useful discussions with R. Benzi, S. 
Succi and R. Verzicco. All numerical simulations were 
performed on the APEmille computer at INFN, Sezione 
di Pisa. 
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